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Q ■ We derive the exact late-time asymptotics for small spherically symmetric solutions 

t^j- ■ of nonlinear wave equations with a potential. The dominant tail is shown to result 

from the competition between linear and nonlinear effects. 



I. INTRODUCTION 



We consider linear and nonlinear wave equations with a potential term 

Uu + XVu = F(u) □ = dl - A , (1) 



>• ' in three spatial dimensions for spherically symmetric initial data 

m : 

^ 1 u(0,r) = /(r), d t u(0,r) = g(r), r := \x\ (2) 

with /, g of compact support. The spherical symmetry of the initial data is preserved in 
evolution so u = u(t,r). We are interested in the asymptotic behaviour of u(t,r) for late 
t^. | times t ^> r. 

Our approach is based on the perturbative calculation which has been developed by the 
last three authors in concrete physical applications [H, 0, 0] and recently put on the rigorous 
ground by the first author in [4| (below referred to as part I). In part I the convergence 
of the perturbation scheme was proved in a weighted space-time L°°-norm which provided 
pointwise estimates on the solution u(t,r) in the whole spacetime. Moreover, upper bounds 
on the errors (remainders of the perturbation series) for every perturbation order were 
obtained. Here, we are going to combine the qualitative global weighted-L°° estimates with 
the quantitative perturbation scheme in order to obtain precise late-time asymptotics of 
solutions. To this end, we first solve the linear perturbation equations analytically up to the 
second (nontrivial) order (in spherical symmetry this can be done explicitly) and show that 
our decay estimate is optimal. Then, we prove that the sum of all higher-order perturbations 
does not modify the dominant asymptotics, hence the second order perturbation gives the 
precise approximation of the tail of the solution u. Along the way, we illustrate our analytical 
results with numerical solutions of the initial value problem ([T]r2]). 

The basis of our analysis is given by the theorem of Strauss and Tsutaya 0], recently 
generalized by NS 0, which states that 

C 

\u(t,x)\ < 7- — , , ,,, : V(t,x) e Mi +3 (3) 

1 V n - (1 +t + |x|)(l + \t- \x\Dv- 1 K ' 1 + K ' 
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with q := min(m — 1, k,p — 1) provided that the potential V and the initial data /, g satisfy 
pointwise bounds 

1/(1)1 £ (TT^ |V/W|£ (TW' |9WI "(TW- m>3,(5) 

with small V , f , fi, g and the analytic nonlinearity satisfying for p > 1 + v2 

\F(u)\ < Fi|u| p , |F(u)-F(v)| < F 2 |w- w|max(|w|, H) p_1 for |u|, |u| < 1. (6) 



This is true for classical solutions [6|, i.e. for (/, g) G C 3 x C 2 , V 6 C 2 and F G C 2 , leading 
to u G C 2 and remains true also for weak solutions, e.g. w G C°, @. 

Here, for simplicity, we consider initial data of compact support so the decay rate q is 
determined solely by the spatial decay rate of the potential k and the leading power of 
the nonlinearity p. Generalization of these results to the initial data with the fall-off (JSJ) is 
straightforward. 

The paper is organized as follows. We first study the purely linear situation with the 
potential term only. Then, we repeat the calculations for the purely nonlinear case without 
the potential. Finally, we combine both results in the general case (JTJ). 



Notation 

We use the symbol (x) := 1 + |x| to denote the spatial weighted-L°° norm 

\\f\\Lg:=\\(r) n f{r)\\u-M (7) 
We also define a space-time weighted-L°° norm 

NU~ := \\(t + r) s (t - r) q ~ s u(t,r)\\ L ^ {R+xR+) . (8) 

We will frequently use the fact that the fmiteness of guarantees the decay of u like 

1/t on the lightcone t ~ r and like l/t q for fixed r as well as l/r q for fixed t. Note that 
functions with compact support in M. + belong to all spaces with any m > 0, what we 
will denote by L^. Analogously Lf^ will stand for functions that belong to Lf q for any q. 

We introduce the following notation for solutions of the wave equations. Let Jy be a 
linear map from the space of initial data to the space of solutions of the wave equation 
CO)-© with F(u) = 0, so that u = Iv(f,g)- For the wave equation with a source term and 
zero initial data 

Du + Vu = F, w(0,r) = 0, <9 4 w(0,r) = 0, (9) 

we denote the solution by u = Ly(F), where Ly is a linear map from the space of source 
functions to the space of solutions to the above problem. Note that, due to linearity, the 
solution u of a wave equation with source F and nonzero initial data /, g is the sum of these 
two contributions 

u = L v (F) + I v (f,g). (10) 
Observe that if we put the potential term on the r.h.s. we obtain 



Uu = -Vu + F, 



(11) 
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which, treated as a wave equation without potential (on the l.h.s.), is formally solved by 

u = -L (Vu) + L (F)+I (f,g). (12) 

Here the solution u appears on both sides what seems to make the formula useless, but it 
will allow us to formulate various iteration schemes, e.g. 

(Vu n )+L (F(u n )) + I (f,g) (13) 

for which we will study convergence in suitable Lf q norms. 

Finally, we define constants which arise from estimates proved in j5| and improved in [9| 



c - :=max (27^)' 5 )' (14) 

C M := 2 + JL + 2 
p — 1 q — 1 

The latter will be referred to as a bound on the allowed strength of the potential. We wish 
to emphasize that this bound, although not optimal, is not arbitrarily small but finite, which 
is crucial in applications (like, for instance, the Regge- Wheeler equation describing waves 
propagating on Schwarzschild geometry). 

We recall some standard definitions of asymptotic calculus. The notation f(t) = O {g(t)) 
for t — > oo means that there exist constants C, T > such that 

\f(t)\<C\g(t)\ (16) 

for all t > T. The notation f(t) = o(h(t)) for t — > oo means that 

lim = 0. (17) 

t^oo h(t) 

We will also use the symbol f(t) = g(t) for an asymptotic approximation, as a shorthand to 
f{t) = g(t)[l + o(l)] as t — > oo. In case when we write fit) = ct~ q and the constant c may 
become zero, this notation should be read as fit) = ct~ q + o{t~ q ). 



II. LINEAR CASE WITH POTENTIAL 

First, we consider the linear wave equation 

Du + \V(r)u = (18) 

with initial data ([2]), where f{r) and g{r) are supported on the interval r G [0,i2]. We 
assume that V(r) = Vo/r k for r ^> 1 and A > is a small parameter, bounded by some 
finite constant Cy > (which will be defined later). Moreover, we assume that the potential 
V and the initial data /, V/, g are (at least) continuous and satisfy 

\\V\\l? = 1 (19) 

and 

/o := H/IU-, ft ■= l|V/|U r+1 , 9o ■= IMU~ +1 (20) 
with f , fi, g < oo for some k > 2. 
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A. Perturbation series 

We define the perturbation series 

AX- (21) 



oo 

U 

n=0 



Inserting ( 1211 into equation ( [TBI we get the following perturbation scheme 

□«o = 0, K^)(0) = (/,s) -H- v = I (f,g) (22) 

□?i n+ i = "^n, (w n+ i,w n+ i)(0) = (0,0) -> = -Lo(^w n ) . (23) 

Due to linearity of ( [18]) it turns out that the partial sums 

n 

u n := A k v k , n > (24) 

fc=0 

satisfy the following iteration scheme 

it_i := (25) 

u n ■= W, 9) - XL (Vu n ^), n > 0. (26) 
Then, from part I, we have the following 

Theorem 1. For f,g and V as above and any k > 2, the sequence u n converges (in norm) 
in L^ k provided that A < C k \. The limit u := \\Ta. n ^ too u n satisfies 

^^ {t + rnt-r}^ V(( - r > < 27) 
where a positive constant C depends only on fo, fi,go, A and k. 
From the proof of theorem [1] it follows that 

\\u n — U n —i 

Klk* = ^ ~ < (C fe)fc )1|/ (/^)||^, (28) 

hence v n G for all n > 0. At the lowest order vq = u$ we have an arbitrarily fast 
decay estimate, Vq G L^^, as follows from Huygens' principle. All higher-order terms 
v n ( n — I; 2, ...) contain contributions from the backscattering off the potential and are only 
in Lf k . Since u G L^ k , we expect that all u n starting from U\ G L^ k predict qualitatively 
correct asymptotic behaviour of u. 

B. Optimal decay estimate 

Theorem 2. Under the above assumptions, for t ^> r + R, we have 

Vl (t,r) = Cl t~\ (29) 
where the constant c\ is given by ( f^Pj) . 
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Proof. For v satisfying (1221) from lemma [T] we have 

Ua{t , r)=M t,r) = h{t - r) - r h{t + r \ (30) 

where /i is given by flA.31) . To solve equation (1231) we use the Duhamel representation for 
the solution of the inhomogeneous equation \3v = N(t, r) with zero initial data 

t t+r-T 

v(t,r) = ^jdr J pN(r,p)dp. (31) 

\t-r-r\ 

This formula can be easily obtained by integrating out the angular variables in the standard 
formula <p — G ret * N where G ret (t,x) = (2Tr)~ 1 Q(t)S(t 2 — \x\ 2 ) is the retarded Green's 
function of the wave operator in 3 + 1 dimensions. It is convenient to express fl3Tl) in terms 
of null coordinates £ = r + p and 77 = r — p 

t+r t—r 

V ^ r) = Tr J * JdV^^-N^v), (32) 

\t-r\ -£ 

where N(£,rj) := N (^p, ^j 1 ) — ^( r ?p)- Using this representation we get from (1231) 

t+r t—r 

v 1 (t,r) = ~ J d£ J dT,£^V(p(t,r,))vo(£,T,). (33) 

\t-r\ -£ 

Since the initial data /, g are supported on [0, R], the function h{x) is supported on [-R, R]. 
Then, for t > r + R equation ( 1331 simplifies to 

+if t+r 

vi(t,r) = ~ Jd v h( v ) J d£V(p(£,T,)). (34) 

—R t-r 

Next, we write the potential in the form V(r) = r~ k [Vo + w(r)] with w (r) — > as r — > 00 at 
any rate (i.e. w(r) = o(l) for r ^> 1). Then 

+if t+r 

Vl (i,r) = ~ y dr/frfo) y ^ [Vb + w(e - 77)] = ^ (t, r) + (t, r) , (35) 



-if, t-r 



where 



+if t+r 

•>k-2 



Mt,r) = -—V J drjh(rj) J di(i- v y k (36) 



-R t-r 
-R t+r 



o 2 p p 

vt(t, r) = -— J dr] h(rj) J d£ (f - v y k w^ - rj). (37) 



-if, t-r 
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Using lemma [2] for t 3> r + R and k > 2 we get 



r + R 
t k+l 



with 



+R 



ci = -2 fc ~ V / %) d/j 



Using lemma 2 again, we get an estimate for v± 



\vi(t,r)\ < 



tk-2 



+R 



t+r 



drj\h(rj)\ sup \w(( - r])\ / d^^-rj)' 

t-r<(<t+r J 
t-r 



< sup 1^(01 

t-r-R<C<t+r+R 



sup moi 

t-r-R<(<t+r+R 



ifc-2 



t+r 



dv\h(v)\ / de(e-»?)' 



t-r 



Cl 



+ 



r + 



where 



ci = -2 



fc-i 



dry. 



Note that the prefactor in (140]) vanishes asymptotically for large i 



lim sup 

t-r-R<C<t+r+fi 



H0\ = o, 



hence, for t > r + i? we have 



r + R 



Cl 







r + R 



Cl 



(38) 



(39) 



(40) 



(41) 



(42) 



(43) 



If the potential behaves like V(r) = V Q r k + W(r) with W(r) — (r k l ) for r > 1, it 
follows from (BSD that 



«i(*,0 = ^ + o 



r + R 
t k+1 
l + r + R 











r + i? 



which gives the more detailed information about the sub-leading term. 
Theorem 3. Under the assumptions of theorem^, for t ^> r + R, we have 

u(t,r) = \ Vl (t,r)[l + 0(\)}, 

hence 



u(t,r) = Ct' 



C = \ Cl + (A 2 ) 



(44) 
□ 

(45) 
(46) 
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Proof. Knowing that the perturbation series converges for some A we can bound the error 
in the n-th order relative to the exact solution by estimating the sum of all higher order 
terms. For the convergent sequence u n we get from the proof of Theorem [1] that 



\U - U n \\ L oo 



< 



(Ck,kX) 



n+1 



1 — Cfc,fcA 

what provides the pointwise bound on the error 



\W,9)\\Lf h , 



i. \ m ^ (Cfc,fcA) n+1 Ck+i ■ (Jo + fi + go) 
u(t, r) — u n (t, r) | < - — — — — • — — 7-r-, — Vi, r > 0. 



l-C k)k X {t + r){t-r) k -* 
For n — 1 with u\ = Vq + Xv± we have 



\u(t,r) - v (t,r) - \vi(t,r)\ < 



(C k , k \) 2 C k+1 -(f + fi + go) 



l-C k>k X (t + r)(t-r)*-i 
A simple inequality (which follows immediately from Bernoulli's inequality) 



: A 1 (t,r) 



< 



implies that 



1-aC 
1 



< 2, VC < l/(2o-), a > 



< 



(47) 



(48) 



(49) 



(50) 



(51) 



for C := r/(l + t) < l/(2q), hence it holds for all t > 2gr. The error term can then be 
estimated by 



A 1 {t,x)<2{C ktk \) 



,2 a 



fc+i • 



(/o + + po) ^ ^ 



(52) 



(l+t) fc 

where we have used twice the inequality fl50|) for t > 2(fc — l)r and A < l/(2C k ,k)- 

From Huygens' principle for ( |22l) with initial data of compact support it follows that 
vo(t, r) = for t > r + R, hence for every r > and sufficiently large t > max[r + i?, 2{k — l)r] 
we have 

fX 2 ' 

\u(t,r) - X Vl (t,r)\ < \u(t,r)-v (t,r) - X Vl (t,r)\ + \v (t,r)\ = O 



(53) 



and 



Cl 



Cl 



- A^ < |w(t,r) - A«i(t,r)| + A ui(t,r) - 



(54) 



where we have used the result of theorem [2], eq. (1381) . Therefore 

u(t,r)^^, C = Xc l + 0(X 2 ). 



(55) 
□ 



This gives the precise quantitative information about the late-time tail of u(t,r) and 
shows that the estimate in theorem [1] is optimal (for t >r) (see Table 1 and Fig. 1 for the 
numerical verification). 



s 







AVb = 10" 3 


Vb 


Theory 


Numerics 


Theory 


Numerics 


k = 3 Exponent 
Amplitude 


3.0 

-3.5449 x 10" 1 


2.9996 
-3.0429 x 10" 1 


3.0 

-3.5449 x 10~ 3 


3.0000 
-3.5394 x 10~ 3 


0.3485 


k = 4 Exponent 
Amplitude 


4.0 

-7.0898 x 10- 1 


4.00001 
-6.6885 x 10- 1 


4.0 

-7.0898 x 10~ 3 


4.00000 
-7.0856 x 10~ 3 


0.2339 


k = 5 Exponent 
Amplitude 


5.0 
-1.4179 


5.00000 
-1.3745 


5.0 

-1.4179 x 10~ 2 


5.00000 
-1.4175 x 10~ 2 


0.1560 



TABLE I: Linear case with a potential term XV (r) = AVb ta ° r k ■ The values of Vb follow from 
the normalization condition (I19p . The results are obtained for the initial data of the form: f(r) = 0, 
g(r) = 4(r 2 — 1) exp (— r 2 ) , which corresponds to h(z) = z 2 exp (— z 2 ) (see lA.3p . The number at the 
Theory- Amplitude entry gives the value of Xc\ with c\ defined in ([39]) . Note that for A Vb = 10 _1 
the values of A are actually greater than the convergence radius of perturbation series obtained in 
theorem [IJ but still, the formula (f55|) seems to work very well. 




FIG. 1: We plot (on log-log scale) the numerical solution u(t, r = 1) of equation (|18p for 
XV (r) = O.ltanh r/r 5 (this corresponds to A = 0.64). The initial data are: f(r) = 0, 
g(r) = 4(r 2 — 1) exp (— r 2 ), which corresponds to h(z) = z 2 exp (— z 2 ) (see IA.3j) . The first three 
terms in the perturbation expansion (|2ip are superimposed. In agreement with theorem 3, the tail 
is perfectly approximated by Xv\ (cf. Table 1). 

III. NONLINEAR CASE WITHOUT A POTENTIAL TERM 

Now, we consider the nonlinear wave equation of the form 

Du = F{u) (56) 

with initial data (/, g) supported on the interval r e [0, R] and satisfying fl20l) with f , f\,go < 
e. The nonlinear term obeys ^(w)] < -^i|^| p for |u| < 1 and \F(u) — F{v)\ < F 2 \u — 
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v | max(|w|, \v\) p 1 . The second condition is satisfied e.g. for F(u) = u p with F 2 = p or for 
F G C 1 such that \F'(u)\ < F 2 \u\ p ~ 1 for \u\ < 1. 

A. Perturbation series 

In order to construct a well-defined perturbation scheme to all orders we have to assume 
additionally that F(u) is analytic at u = and its Taylor series starts at power p > 3 

oo 

F(u) = u P Y, b nU n , b ^0 (57) 

n=0 

Then, for small initial data 

(u,u)(0) = (ef,eg) (58) 
we introduce the perturbation series for the solution of (156]) 

oo 

u = £Vi; n - (59) 

71=1 

Inserting this series into 056|) and collecting terms according to powers of e we obtain the 
following perturbation scheme 

□«i = 0, («i,*i)(0) = (/,^) - «i = io(/,y) (60) 

□un+i = F n («!,... ,v n ), (w n+ i,?) n+ i)(0) = (0,0) -> = L (F„(wi, ...,v n )), (61) 

for n > 1, where _F n result from collecting the nonlinear terms with the same powers of e 

n , 1 n , n 

F n ( Vl ,...,v n )J2^>x k ■■■vl k , (62) 

k 

where a^' m G N satisfy J^ =1 ^Q!fc' m = n + 1 and Xlm=i a fc' m — P f° r ever Y n i k- The 
coefficients are functions of b m only (see 0] for the explicit formula). 

We call this expansion the "zero background" case because the zero-order term v is 
absent. If a Vq term was present in the series above (i.e. the summation started at n — 0), 
we would have an additional equation Dvq = F(vq) which is genuinely nonlinear (in contrast 
to the above system of linear wave equations with source terms) . Its solution vo represents 
a "background" around which the perturbations v n are calculated. 

From part I, we have the following 

Theorem 4. With f, g and F(u) as above for any p > 3 and sufficiently small e the series 
defined in / f5Pj) -/ f^7]) converges (in norm) in Lf p _ l to the solution of equation l[5b]) with initial 
data (ESP- 

Since the introduction of the auxiliary parameter e in the perturbation series expansion 
serves only to generate a system of linear equations equivalent to the original nonlinear 
equation, we can eventually remove the parameter e and assume that the initial data are 
such that /o,/i,<7o are sufficiently small. Then, solving the system (!60l) - (l6Tl) and summing 
up the convergent series Y^=i v n u we obtain a solution of the nonlinear wave equation (!56l) . 
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B. Optimal decay estimate 

The perturbation scheme fl60l) - fl6Tl) can be written as 

vi=I (f,g) (63) 

v 2 = v 3 = ... = v p _i = (64) 

v p = L (F p _i(vx, Vp-i)) = b L ((v 1 ) p ) (65) 

v n +i = L (F n (vi, v n )), n>p. (66) 

We have V\ = I (f, g) G and v n G L% x for n > 2. 

Theorem 5. Under the above assumptions, for t ^> r + R, we have 

v p (t,r)^d p r ( P- 1 \ (67) 

where the constant d p is given by [72$ . 

Proof. In analogy with equations (I30H32P we have from lemma [1] 

, , hit - r) - hit + r) 
v 1 (t,r) = -± '-^ (68) 

and from (|65p 

t+r t—r 



v p (t,r) 



1 J di J dviZ-rifoiv^Oy. (69) 



\t-r\ -£ 

As before, interchanging the order of integration we get for t > r + R 

+R t+r 

-3 



v p (t, r) = ^— 1° J dr, (h( V )Y J dat- V)~ P+1 - (70) 



-R t-r 

Using lemma [2] we get for t ^> r + R and p > 3 



t p-i y t p 

where 

+R 



d p = 2^X J d V (h( V )y. (72) 

—R 

□ 

Now, we will show that v p dominates the perturbation series for large times and small 
e and has the same decay rate as the full solution u of the nonlinear wave equation (see 
Table 2 and Fig. 2 for the numerical verification). 

Theorem 6. Under the assumptions of theorem^ for small e and t ^> r + R, we have 

u(t,r)=e p Vp(t,r)[l + 0(e)], (73) 

hence 

u(t, r) ^ Dt~ p+1 , D = d v e v + O (e p+1 ) . (74) 
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Proof. We need to show that eIo(f, g) and e n+1 L (F n (vi, ...,v n )) for n > p are small relative 
to e p d p (x)t~( p - 1 ) . As before, for Vi = I (f,g) E L^ 3 ^ Huygens' principle and compact 
support of the initial data imply that Vi(t,r) = for sufficiently large t (and fixed r). 
From the convergence proof for the perturbation series we know that there exist constants 
M,p > such that ||tVi||Lj° _ < Mp n for all n > 1. Hence, for sufficiently small e < 1/p we 
can estimate the remainder of the perturbation series 



m=p+l 



e m v r 



i.p-i 



<M y, 

m=p+l 



e p 



3+1^+1 

1 — ep 



This implies that for t ^> r 



e m v m {t,r) 

m=p+l 



< 



Ce p+l 



(t + r)(t — r)* 



O 



(75) 



(76) 



so 



\u(t,r) - e p v p (t,r)\ < \ev 1 (t,r) \ + 



m=p+l 



o 



r P+l 



From theorem for t ^> r + R we have = dpt'^-V + (r p ) , hence 











which finally gives 



u(t, r) Dr p+ \ D = d p e p + O (e p+1 ) 



(77) 



(78) 

(79) 
□ 





e = 1 


e = 10 _1 


Theory 


Numerics 


Theory 


Numerics 


p = 3 Exponent 
Amplitude 


2.0 
0.1421 


2.0009 
0.1265 


2.0 

0.1421 x 10" 4 


2.0008 
0.1427 x 10" 4 


p = 4 Exponent 
Amplitude 


3.0 

9.0873 x 10~ 2 


3.0013 
8.4433 x 10" 2 


3.0 

9.0873 x 10~ 6 


3.0012 
9.1631 x 10" 6 


p = 5 Exponent 
Amplitude 


4.0 

5.9925 x 10" 2 


4.0015 
6.1192 x 10" 2 


4.0 

5.9925 x 10" 7 


4.0015 
6.0597 x 10" 7 



TABLE II: Nonlinear case: F(u) = u p , without a potential term. The initial data are the same as 
for Table 1 and Fig. 1. The number at the Theory- Amplitude entry gives the value of e p d p , with 
dp defined in ([72]) (with 6q = !)• 



IV. NONLINEAR CASE WITH A POTENTIAL TERM 



Finally, let us consider the full nonlinear wave equation (TjQ) with a potential with initial 
data (/, g) supported on the interval r e [0, R] and satisfying (120]) with fo,fi,go < e - The 
nonlinear term F(u) is the same as in the previous section. 
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1 

0.01 



1e-04 



p- 1 e-06 
j[ 

^ 1 e-08 
1e-10 
1e-12 
1e-14 

0.1 1 10 100 1000 

t 

FIG. 2: We plot (on log-log scale) the numerical solution u(t, r = 1) of equation (|56p with -F(ii) = 
u 3 . The initial data are the same as in Fig. 1 and e = 0.1. The first three terms in the perturbation 
expansion (|57h are superimposed. In agreement with theorem 5, the tail is perfectly approximated 
by e 3 v 3 (cf. Table 2). 

A. Perturbation series 

Defining the perturbation expansion for the nonlinear wave equation with a potential (CEJ) 

oo 

u = £)eX, (80) 

n=l 

we encounter the problem of two scales which are given by parameters A (measuring the 
strength of the potential) and e (measuring the strength of the initial data). Since these 
parameters play only an auxiliary role in generating the perturbation scheme, we make a 
convenient choice and assign to A a scale of some power of e, say A = \e a with a G N + . 
Then, the power series ( 1801) inserted into the wave equation §S§ gives 

w_ n :=0, n>0, (81) 
v 1 :=I (f i g), (82) 

v n+1 := -\L (Vv n+1 - a ) + L (F n (vi, ...,v n )), n>l. (83) 

In the following we choose a := p — 1 because then the lowest-order nontrivial term, v p (all 
lower-order terms with 1 < n < p vanish), contains contributions both from V and F and 
gives a good approximation to u, as will be shown below. 
In this case, from part I, we also have a convergence result 

Theorem 7. With f, g, V and F(u) as above for any k > 2, p > 3, A < C~l and sufficiently 
small e the series defined in (E5P-(E^ converges (in norm) in L™ q for q = min(p — 1, k) to 
the solution of the equation (fip with initial data |53J). 




13 



B. Optimal decay estimate 



For a = p — 1 the system fl8Tj) - fl83]) takes the form 

v- n := 0, n > (84) 

vi = I (f,g) (85) 

V2 = v 3 = ... = v p ^x = (86) 

v p = -XL (Vvi) + L (F p _i(vx, t'p-i)) = -XL Q (V Vl ) + & A)(0i) p ) (87) 

v n+1 = -\L (Vv n - p+2 ) + L (F n (vx, v n )), n> p. (88) 
Theorem 8. Under the above assumptions, for t 3> r + R, we have 

v p (t,r) = dpt~ q , g :=min(A;,p- 1), (89) 



where the constant e p is defined in (Wfy. 

Proof. v p defined in equation ( 1H71) is a sum of two contributions, from the potential and from 
the nonlinear term, 

Vp(t, r )^v^(t, r ) + v; on (t,r), (90) 

where 

v%°%r) := -\L (V Vl ), v n ™(t,r) := 6 ^oK)- (91) 
From theorem [2] we have 

<%r) = % + o(±-) (92) 



p v ' ' t k \t k 

with 



+R 



Cp = -2 k - l \V I h(ri)dri, (93) 

and from theorem [5] we have 



-R 



t p-i y t p 

with 

+R 



d p = v-x j d v (h( v )y. (95) 



-R 



Depending on whether k<p— loik>p— 1, the linear v pot (t, r) or the nonlinear Vp° n (t, r) 
contribution to the tail is dominant, respectively. In the special case k = p — 1 we have 



Thus, the constant in (I8"9~j) is given by 



c p if k < p — 1 

c p + dp if k =p - 1, (97) 
dp if k > p — 1. 

□ 
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Now, we will show that v p dominates the perturbation series for large times and small e 
and has the same decay rate as the full solution u of the nonlinear wave equation with the 
potential (see Fig. 3 for the numerical verification). 



1e-14 



1e-16 



linear 
non-linear 
full solution 




FIG. 3: We plot (on log-log scale) the numerical solution u(t, r = 1) of equation (TT]) with F(u) = u 3 . 
The potential AFand the initial data (ef, eg) are the same as in Fig. 1 with A = 0.64 and e = 0.001. 
Superimposed are solutions with the nonlinearity or the potential switched off. The crossover from 
the linear tail ~ i -5 (for intermediate times) to the final nonlinear tail ~ t~ 2 is clearly seen. 



Theorem 9. Under the assumptions of theorem^ for small e and t ^> r + R, we have 

u(t,r)=s p v p (t,r)[l + 0(e)], (98) 

hence 

u(t,r)^Et- q , q := min(fc,p- 1), E = e p e p + O (e p+1 ) . (99) 

Proof. We can repeat the reasoning from the proof of theorem [6] where we used the fact that 
the perturbation series ^ n=1 ^ n i'n is convergent. Here, theorem [7] guarantees convergence in 
L^° q with q := min(/c,p — 1). Analogously, we obtain for t 3> r 



< 



Ce p+1 



(t + r) (t - r)*- 1 



n=p+l 

so (again, Vx(t,r) vanishes for t ^> r by Huygens' principle 
\u(t, r) - e p v p (t,r)\ < |ei; 1 (t,r)| + 



O 



ifl 



(100) 



£my - 

m=p+l 



o 



rP+1 



(101) 



From theorem [8] we have v p = e p t q + o(t q ) for t ^> r + R, so we get 



u(t, r) — 



e p s p 



t" 



O 



^) +0 (£! 

ti \ti 



(102) 
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which gives 

u(t, r) = Et~ q , E = e p e p + O (e p+1 ) . (103) 

□ 
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APPENDIX: LEMMAS 

Lemma 1. The solution of the free wave equation 

Du = (A.l) 

with spherically symmetric initial data u(0,r) = f(r), d t u(0,r) = g(r) has the form 

, , hit - r) - hit + r) . A , 

u (t, r) = -± > ^ i, A.2 

r 

where 

h(r) = - r -f(r) + ^J°° r'g(r')dr', (A.3) 

which is defined for all r G K. by the extension /(— r) := fir), g(—r) := g{r). When f and 
g have compact support then h has also compact support on R. 

Proof. In spherical symmetry the wave equation (lA.lj) can be written as 

O^iru) = (A.4) 

where £ = t + r and rj = t — r. Its most general solution has the form 

ru(t, r) = h-{rf) + h+(£) = h-(t-r) + h+(t + r). (A.5) 

We require that u(t,r) be finite at r = what implies 

= h-{t) + h+{t) h{t) := h-{t) = -h+{t). (A.6) 

From the initial conditions we get 

/(r) = u(0 , r)= M-r)-fc(r) (A?) 

r 

^(r) = 8 t u(0, r) = h '(- r )~ h '( r \ (A.8) 



We can write 



Kr) = \[h{r) + h(-r)} + l -[h{r) - h(-r)} = ±S(r) + ^A(r), (A.9) 
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where we have introduced a symmetric function S(r) := h(r) + h(—r) and an anti-symmetric 
function A(r) := h(r) — h(—r). The solutions for A(r) and S'(r) = h'(r) — h'(—r) can be 
immediately read off from the initial conditions (1A.71 IA.8j) : 



A(r) = —rf(r) 
S'(r) = —rg(r). 



(A.10) 
(A.ll) 



We see that the extension of / and g on all r £ R defined by f(—r) '■— f(r), g{—r) := ^(r) 
is the consistency condition for eqs. (1A.101 [A. lip . Integrating (lA.lip we get 



S(r) - S(0) 



r'g(r')dr'. 



We use the freedom of choosing the integration constant and set 



(A.12) 



(A.13) 



what gives (1A.3I) . With this choice we obtain h(r) compactly supported on R if f(r) and 
g(r) are compactly supported. To see this, assume f(x) — g(x) = for \x\ > R and consider 
r > R. The function h(r) is obviously zero from (1A.3I) . For negative arguments 



S(0) ■= I r'g(r')dr' 



h(—r) = - f(—r) H — / r'g(r')dr' = - r'g(r')dr' = 



R 
=0 

because the integrand r'g{r') is an odd function. Thus, supp h C [-R, +R] 
Lemma 2. Let a > 1. T/ien 

+ R t+T 

h{rf) drj 



-R 



t-r 



-/ h(rj)drj + 

1 J-R 



r{r + i?) 



(A. 14) 



□ 



(A.15) 



for t > lair + R) and all r > 0. 

Proof. Consider first the inner integral for 7] £ [— i?, -R] 

"* +r di f +r dy 1 



+r 



dy 



2r 1 



+ r 



We have 



+ r 



2r 1 



(A.16) 



|<f(t,r,77)| < 
and the integrand J can be estimated by 



r/.iy 



(A.17) 



J := 



- 1 



< 



(A. 18) 
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what can be shown as follows. Having in mind that — (r + R) <y + r]<r + Rwe find for 
y — r] < 

j = - L___i<_ L——-!^^ (A.19) 



and for y — r] > 



J = 1 



(l + ita)° 
1 



< l 



J 2 



(A.20) 



By simple algebra one can easily show that J2 < J\ what gives (1A.18I) . Further, by a version 
of Bernoulli's inequality, 



r+R 



1 a 

J < — ~ 1 * 



Then 



a 



r+R 



W,r,r,)\ < Y^M 2r ^ 4a^±^ 



for t > 2a (r + R). Finally, 



+R 



h(rj) drj 



t+r 



t-r 



+R 



h(rj) dr) 



R 



2r 5(t,r,rj) 



2r 



+R 



h(rj) dr] + 



r(r + R) 
t a+1 



(A.21) 



(A.22) 



(A.23) 



□ 
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